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We investigate the possibility of using dielectric microdisk resonators for the optical detection 
of single atoms trapped and cooled in magnetic microtraps near the surface of a substrate. The 
bound and evanescent fields of optical whispering gallery modes are calculated and the coupling 
to straight waveguides is investigated using finite-difference time domain solutions of Maxwell's 
equations. Results are compared with semi-analytical solutions based on coupled mode theory. 
We discuss atom detection efficiencies and the feasibility of non-destructive measurements in such 
a system depending on key parameters such as disk size, disk-waveguide coupling, and scattering 
losses. 
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I. INTRODUCTION 

Micro resonators are a promising system for research 
in a wide range of fields. Their spectral properties can be 
exploited in various applications ranging from telecom- 
munication to biological/chemical sensors 0] as well 
as in fundamental research. Of specific interest is the po- 
tential contribution of such devices to the emerging field 
of quantum technology (QT), which may again serve as 
an enabling technology for both fundamental science and 
applicative research. 

In the realm of QT, single atoms and ions coupled to 
micro resonators are one of the most promising systems. 
There, a high degree of control over light-atom interac- 
tion can be achieved, which may enable new insight and 
capabilities in contexts such as cavity QED single 
photon sources memory, and purifiers for quantum 
communication j5|, manipulation of internal and exter- 
nal degrees of freedom for matter waves in systems like 
interferometric sensors 0, clocks, and perhaps also the 
elusive quantum computer 0, Q • In this paper we will 
analyze a specific micro resonator in the form of a mi- 
crodisk, and furthermore we will focus on its ability to 
detect single atoms. 

Recently reported progress in the manufacturing of 
high Q dielectric microdisk structures '3| motivates the 
development of compact and integrable devices. Of spe- 
cific interest is the wafer-based manufacturing of res- 
onators where a good control of the physical characteris- 
tics can be achieved during fabrication enabling, for ex- 
ample, extremely small mode volumes as well as accurate 
alignment with other elements such as microtraps. Fur- 
thermore, a wafer based device may allow more complex 
functions such as tunability to be integrated. 

Experiments have also shown that it is possible to trap, 
guide and manipulate cold, neutral atoms in miniaturized 
magnetic traps above a substrate using either microscopic 
patterns of permanent magnetization in a film or mi- 
crofabricated wire structures carrying current or charge 
Such surfaces have received the name atom chip. 



Recently it was proposed to integrate micro-optics into 
an atom chip for ato m-lig ht interaction. In particular, 
the use of Fabry-Perot [Tj or photonic bandgap 12] cav- 
ities has been discussed. As a complementary and com- 
parative analysis we investigate in this work the external 
fields of optical whispering gallery modes which are sup- 
ported by a toroidal microcavity and which may be used 
for the above purpose. Such a microdisk can be consid- 
ered as a two-dimensional version of the much studied 
microsphere [l3llT3.ll^ . 

This work is organized as follows: In Sec^jwe describe 
the system under consideration. Section IIIII deals with 
the different numerical and analytical methods which we 
apply to examine the disk resonator structures. We an- 
alyze feasible realizations, taking into account imperfec- 
tions of various origins. We then introduce atoms to the 
system in the framework of the Jaynes-Cummings model 
(Sec. II V|) . Subsequently, the calculated fields are used 
to estimate the efficiency of our scheme for single atom 
detection (Sec. |V||. In Sec. EH we discuss experimental 
considerations such as tunablity. Finally, we conclude in 
SeclVTH 



II. SYSTEM DESIGN 

The basic system under discussion is that of an atom 
chip consisting of a magnetic or dipole microtrap for cold 
atoms, and a photonic part which is used for optical atom 
detection. As shown in Fig. ^ in this work we focus on a 
photonic system consisting of a circular plate (disk) and 
a linear waveguide. 

The waveguide couples light into and out of the disk. 
In a realistic setup, both waveguide ends will be attached 
to optical fibers. In order to optimize power transfer be- 
tween the waveguide and a single-mode fiber, the mode 
overlap at the interface has to be maximized. This re- 
quires waveguide dimensions of 9x9 to 12x12 microns in 
the refractive index range of 1.454 - 2.17 for wavelengths 
around 780nm (the wavelength of Rb atoms usually used 
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FIG. 1: The structure under consideration. The evanescent 
wave from the slab waveguide (1) is coupled into the disk 
(2) and back through a small gap between them. The adia- 
batic waveguide tapers (3) serve for coupling light from optical 
fibers (not shown) into and out of the waveguide. Cold atoms 
(4) can be positioned on the disk side. 

in this kind of atom chip experiments). In this case, 
best coupling efHciencies of the order of 96-98% can be 
achieved. For best mode matching with the microdisk 
modes, the waveguide width has to be reduced to about 
0.3 to 1.2/im using adiabatic tapers, as shown in Fig.^ 

The disk itself supports a large range of resonant 
modes. Here, we are mainly interested in the low-loss 
modes traveling along the disk edge in the form of whis- 
pering gallery modes (WGM). While most of the mode 
energy is confined within the disk, a small part of the 
mode exists outside the disk as an evanescent field, and 
it is here that an atom can interact with the light. This 
atom-light coupling changes the optical properties of 
the disk mode, which subsequently changes the inten- 
sity and/or phase of the light at the output of the linear 
waveguide. These changes can then be measured to infer 
the presence of the atom. 

In the horizontal direction the linear waveguide and 
the disk are bordered by air or vacuum with refractive 
index rid = 1. In our calculations we assumed structures 
made of fused silica with ric = 1.454 at a wavelength 
of 780nm. In the vertical direction the structure may 
be more complicated with several layers in order to give 
good mode confinement. 

A number of numerical and analytical tools have been 
utilized to investigate the optical properties of such sys - 
tems, for example coupled mode theory (GMT) |ig.[l7l |. 
scattering matrix theory finite element [T9 |. an d fi- 
nite difference time domain (FDTD) methods [23|. In 
this work we will use FDTD simulations which provide 
rigorous numerical results but which are rather time- 
consuming and therefore not adapt to scan large parame- 
ter ranges. We will thus resort to a semi-analytical CMT 
as a fast tool for a detailed design parameter analysis. 



These two methods complement one another and give a 
powerful tool for the investigation of micro resonators. In 
particular, we are interested in the system characteristics 
dependent on disk diameter, waveguide width, gap width 
between the linear waveguide and the disk resonator, and 
the surface quality of the disk. 



III. OPTICAL PROPERTIES OF 
WAVEGUIDE-COUPLED WHISPERING 
GALLERY MODES 

A. Finite difference time domain calculations 

The first method we use to investigate the resonance 
behavior of the microdisk is by FDTD simulations. Here, 
Maxwell's equations are discretized in space and time, 
and the time evolution of the electromagnetic fields is 
numerically calculated on this grid. To obtain high accu- 
racy it is necessary to make the cells much smaller than 
the optical wavelength, which leads to long calculation 
times, in particular for a three-dimensional (3D) model. 

In two dimensions it is possible to perform direct 
FDTD simulations for relatively large disks of diameter 
> 50fim. In 3D, simulations are only feasible for small 
disk diameters. However, comparisons of 2D and 3D sim- 
ulations for small disk diameters have shown reasonable 
agreement. In the rest of this paper, we will thus re- 
strict ourselves to simplified 2D calculations, where the 
disk and waveguide modes are calculated for a geometry 
infinitely extended in the vertical (y) direction. The 3D 
modes are assumed to be simple slices of thickness dy 
of these infinite modes. For the calculations presented 
in this work we always assume dy — S/im. We will also 
limit the calculations to electric fields polarized along the 
y direction, i.e., to TE modes only. 

In this case. Maxwell's equations for the electric and 
magnetic fields, Ey, H^, and reduce to 

dH, 1 dEy 

1^ = ^ 

dH^^ ldEy^ ^2) 

dt dx ' 

dEy _ 1 OH, dH^ 

where /^o and e are the magnetic and electric perme- 
ability, respectively. The FDTD simulations solve Eqs. 
(^3)-© on a spatial grid for a given initial field distri- 
bution. We used the so-called unsplit perfectly matched 
layer (UPLM) boundary conditions |22 and a uniform 
discretization with a step size range of 0.02/im to 0.08/im 
in the x and z directions and with a time step range of 
4.45 X lO^i^s to 1.78 x lO^^'^s. 

For our calculations the initial condition was such that 
the disk was empty and that light was pumped into the 
lowest transverse mode of the linear waveguide from one 
end. The incoming light was either a continuous wave 
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or a BOfs Gaussian pulse. In the former case, we are in- 
terested in the steady state field distribution which, for 
example, allows us to observe the resonant disk mode 
and the evanescent field. Pulsed input allows the calcu- 
lation of the output power at the other end of the linear 
waveguide as a function of frequency, i.e., the transmis- 
sion spectrum, by applying a discrete Fourier transform 
on the output field, calculating the Poynting power den- 
sity and integrating over the waveguide cross-section. 



B. Coupled mode theory 

We complement the numerical FDTD results with 
a semi-analytical GMT. For this, it is assumed that 
the linear waveguide supports only a single transverse 
mode, while the disk supports two degenerate, counter- 
propagating WGM. The pumped waveguide mode only 
couples to the forward propagating disk mode, but scat- 
tering due to side-wall roughness may couple light into 
the second mode. Both WGM are coupled to the waveg- 
uide mode via their evanescent fields at point 1 in Fig. 
n For the calculations we closely follow the work by 
Rowland and Love |0| . 

First, the WGM are obtained as solutions of Eqs. QJ- 
(Pl in cylindrical coordinates. This gives mode functions 
of the form 



Jl(kncr)_ ±il<t> 



for r < R 



where Ji are Bessel functions and H^^'' are Hankel func- 
tions of the first kind, and R is the disk radius. The 
eigenvalue equation for these modes is given by 

Ji+i{kncR) _ Hj:l\{knciR) 



(4) 



JiikucR) 



rid 



H\^\kn,iR) 



As all WGM are lossy, the eigenvalues are complex 



(5) 



(6) 



and the intrinsic quality factor of the WGM is given by 



iWGM 



'2ki 



(7) 



Similarly, the mode functions Enn of the linear waveguide 
are calculated for the same wavenumber k. 

The second step of the GMT is to write the total light 
field as a superposition of disk and waveguide mode 

E{x, z) = ai{z)Eun{x, z) -f a2[z)EwGM{x, z). (8) 

The coupled mode equations read 

dai 



dz 
da2 
dz 



-i(3unai + iC{z)a2, (9) 
-iPwGMiz)a2 + iC{z)ai, (10) 



where Pun and f3wGM are the propagation constants of 
the waveguide and the WGM, respectively, and C{z) is 
the position-dependent coupling coefficient obtained by 
calculating the mode overlap of waveguide and disk. For 
details of this calculations see Ref. Ilj. 

Finally, integrating Eqs. © and H1Q|I over z in the 
region of significant coupling around point 1 in Fig. ^ 
yields the coupler transmission matrix T which relates 
the waveguide and disk output fields to the inputs. 



where 



ai 
a2 



T = 



T 



ai 

02 



til ti2 
t21 t22 



(11) 



(12) 



The cavity decay rate (half width at half maximum) 
kt of the WGM due to the coupling to the waveguide is 



KT = |ii2|V(2T.) 



(13) 



where Tr = 2-Kl/[ck) is the round trip time. The corre- 
sponding quality factor is 



Icoup = cA:/(2kt) = 27rZ/|ti2p. 



C. System losses 



(14) 



Apart from the intrinsic WGM losses (0 and the cou- 
pling losses into the waveguide (|14|) at least two other 
loss mechanisms have to be taken into account, intrinsic 
material losses and surface scattering losses. 

The main material loss mechanisms are bulk Rayleigh 
scattering and ultraviolet and infrared absorption. The 
corresponding quality factor Qmat can be expressed in 
the form I2II2J 



aA ' 



(15) 



where a is the loss coefficient and A = 27r/fc is the vac- 
uum wavelength. Material losses for fused silica in the 
wavelength range near 780nm are of the order of 5dB /km, 
which gives rise to Qmat 10"^°. 

Greater uncertainty is associated with losses due to 
surface scattering and absorption due to surface rough- 
ness or the presence of an absorbing impurity on the 
surface of the disk. For a given size of surface inhomo- 
geneities (roughness) and correlation length, the surface 
scattering quality factor Qsurf has to take into account 
not only the direct scattering of light out of the disk 
("leakage") but also scattering into other modes with 
high rates of leakage. Various expressions have been 
used to describe quality limits due to surface scattering 
[2^ I2M 01 • In this work, we apply the expression based 
on the model of Rayleigh scattering by molecule-sized 
surface clusters l23l 



Qsurf 



2ic7r2cr2 ' 



(16) 
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where D is the disk diameter, <j is the root-mean-square 
of the surface roughness and the surface correlation 
length. As was reported in the numerical values for 
a and Lc can be less than Inm and 5nm, respectively. In 
our calculations we used cr=lnm, Lc=5nm and cr=2nm, 
Lc=10nm. 

The overall cavity quality factor taking all the loss 
mechanisms discussed above into account is then given 
by 

Q ^ Qcoup ^ QwGM ^ Qmat ^ Q surf ^'^'^) 

This is related to the cavity line width (HWHM) k by 

Q = cfc/(2K) (18) 
and to the cavity finesse F by 

F^Q^^Q/l (19) 
cfc 

where FSR m ck/l denotes the free spectral range {I is the 
mode index of the WGM) . For later use we also introduce 
i^ioss = K — kt which is the cavity loss rate due to losses 
into all channels apart from the linear waveguide. 

D. Results 

After describing the building blocks of our calculation, 
let us now discuss some numerical results for the optical 
properties of our system. This will serve two purposes: 
first, to compare our different calculation methods with 
each other and with available experimental data in or- 
der to show the accuracy of our results; second, we need 
to apply our calculation to the experimental parameters 
that are of interest in our case, in order to establish a 
base for the atom-light interaction that will be discussed 
in the following sections. 

We have compared the spatial profile of the WGM 
obtained from FDTD simulation with the one resulting 
from a CMT calculation, and found good agreement. We 
have also looked at the resonance spectrum of the disk 
for different parameters, and extracted values for the free 
spectral range (FSR) and the quality factor Q. As pre- 
sented below, here too the agreement was good. 

In Fig. 121 we plot the FSR for disk diameters D in 
the range of 5 - 50/im, refractive index values 1.454 - 3.2 
and wavelengths near 0.78yLtm and 1.55/Ltm. The latter is 
chosen to compare with previously published results '2(^ . 
FDTD simulations and analytically calculated WGM are 
in excellent agreement with each other and with the ex- 
perimental data. As expected, the FSR is approximately 
inversely proportional to the disk diameter. 

Figure 13 shows the total cavity quality factor Q, Eq. 
(|17|l . dependent on the outer diameter D of the disk and 
the size of the air gap between disk and linear waveguide. 
Q was calculated using FDTD and CMT models, which 
again show excellent agreement. We note that very high 
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FIG. 2; Free spectral range versus disk diameter. The 
lines present results of analytical calculations at wavelengths 
780nm (solid line) and 1550nm (dashed) in fused silica and at 
1550nm for — 3.2 (dotted). Corresponding FDTD results 
are indicated by (*), and experimental data [2(1 1 by (o). 
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FIG. 3: Quality factor Q versus disk diameter for various 
gap sizes and materials. Solid curves (from bottom to top); 
CMT results for gap sizes O.l/xm, 0.2/im, 0.5/im, and for the 
uncoupled disk for A = 780nm, ric = 1.454, a — Inm, Lc — 
5nm, (•) represent FDTD simulations. Dashed curve: ric = 
2.17, gap size 0.2^im. CMT results for A = 1550nm, a = 
2nm, Lc = lOnm: uncoupled disk with tIc = 1.444 (dash- 
dotted curve), He — 3.2 and gap size O.l^m (dotted curve). 
Corresponding experimental data is given by (o) Q and (o) 
[2o|. respectively. 



quality factors up to about 10^ can be achieved with the 
current system for disk diameters of several tens of mi- 
crons. The upper limit for Q shown in Fig.|31is the value 
obtained for a very large air gap, where coupling losses 
Kt are negligible and the cavity quality is limited by the 
intrinsic and material quality factors QwgMi Qmat, and 

Qsurf • 

Figure O also shows data for two different wavelengths 
and varying refractive index. Note that Q generally in- 
creases with increasing refractive index due to the better 
confinement and therefore weaker coupling to the waveg- 
uide. Similarly, Q decreases with increasing wavelength 
due to stronger waveguide coupling. Nevertheless, for 
each case there are unique material losses and fabrication 
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TABLE I: Optical properties of selected WGM. Qi (Q2) is 
the quality factor Q for a gap size of 0.3/im (0.6^m), / and 
q are the longitudinal and radial mode index, respectively, go 
is the single-photon Rabi frequency for an atom at the disk 
boundary. Surface parameters are a — 2nm, Lc — lOnm. 
Results are obtained using GMT. 



uncertainties which may change this general behavior. 

To calculate the Q factor for different disk sizes, the 
index I of the WGM must be changed accordingly to 
keep the resonant wavelength near 780nm relevant for Rb 
atoms. The necessary wavelength for optimal mode reso- 
nance can be achieved by choosing a precise disk size. We 
will discuss tuning of the micro resonator later in Sec. lVll 
Note also that, while the intrinsic quality factor Qwgai 
is always lower for higher radial modes, the overall Q 
can in fact be higher due to weaker coupling (smaller 
losses) to the linear waveguide, which is the dominating 
loss mechanism for most parameter regimes we are inter- 
ested in. Table shows the optical properties of a few 
selected disk- waveguide geometries. 



IV. THE JAYNES-CUMMINGS MODEL 

In this section we introduce a model for the coupling 
of a two-level atom to the light field outside of the mi- 
crodisk. We consider two coherent modes with com- 
plex amplitudes a+ and a- traveling in opposing direc- 
tions, and assume atom-light coupling to these modes via 
the respective single-photon Rabi frequencies g+ and g- . 
The coupling of the two modes due to imperfections in 
the disk is accounted for by the introduction of the com- 
plex coefficient e, and the model includes the possibility 
of pumping from either direction with rates ry+ and ?/_ , 
where ri± = ti\Ain^±l \/T\.. Here Ain.± is the amplitude 
of the pump field in each direction within the waveguide 
such that |Ami,±P is the photon fiux in units of photons 
per second, and ti\ and have been defined in Scc. lIIIBl 

Ignoring the external motion of the atom, the Hamil- 
tonian of this system can be written as (?i = 1) 



-i(5_aLcroi - .9* o-iofl-) - iri-{a^ 
— i(ea^a_ — e*aLa+) , 



(20) 



where and Ac are the atomic and resonator detunings, 
respectively, ctio and ctqi the atomic raising and lowering 



operators, and aj_ and a± the mode creation and anni- 
hilation operators. Here the energy of the lower atomic 
state has been set to zero. The equation of motion for 
the density operator of the system can be written as 



-i[H,p] + Lp 



(21) 



where L is the usual linear operator describing cavity and 
atomic decay with rates k and F respectively. 

Assuming a factorized density operator p, we find the 
equations of motion for the elements of the atomic den- 
sity operator 



dt 
di 



Pw = 



Pii 



(-r + iAa)piQ 

+{gla++g*_a.){poo- Pii), (22) 

-2Tpii + {g*^a+ + g*_a^)pm 

+ {g+a\+ g^a*_)piQ, (23) 



and the coherent state amplitudes a± obey the equations 
of motion 



dt 
d 



= (iAc - K)a+ - 5+/910 + 77+ - (24) 
— (iAc — K)a- — g-pio + J]- + e*a+. (25) 



In this work we are only interested in the stationary 
solution of these equations of motion. To this end, we 
first solve the linear set of equations H24|l and H25|) with 
respect to a± . The resulting expressions for a± are linear 
in pio and can be inserted into (|22|l . From this and using 
Pii + Poo = 1, we obtain pio as a function of pu. That 
and the corresponding results for a± can be inserted into 
(|23|l to give a real- valued nonlinear equation in pn which 
can be solved by standard numerical techniques. The 
output field of the linear waveguide is the superposition 
of the input field transmitted through the waveguide-disk 
coupler and the light coupled out of the disk. 



Aout,± — tllAin ± 



t 



12 



■.a±. 



The Rabi frequency g± is given by 



.9± 



^WGM 



(Xa) 



3Fc^ 



uj'^dy J rn{r)'^\EwGM{rWdr 
He (rici) for r < i? (r > i?), d, 



(26) 



1/2 



, (27) 

where n{r) — ric (rici) for r < R [r > R), dy is the 
disk height, uj is the atomic transition frequency, and 
Xa — (fa, 0a) is the atomic position in the evanescent 
field of the disk modes. Figure 01 shows the Rabi fre- 
quency depending on the distance of the atom from the 
disk surface for several disk diameters and mode num- 
bers. The maximum values for g± at the disk surface for 
several selected modes are also given in Table U 

We note that for efficient atom-mode coupling the 
atom should be placed within 50-lOOnm from the disk. 
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FIG. 4: Rabi frequency versus atom-disk distance. Solid 
curve: D — 30/im, mode indices {l,q) — (167,1); dot- 
ted: D = SO/xm, {l,q) = (159,2); dashed: D = 15/im, 
{l,q) = (81,1); dash-dotted: D = 45/im, {l,q) = (253,1). 



Such close proximity should be feasible without the atom 
being absorbed by the surface, as the van der Waals forces 
may be counter balanced by blue detuning of the res- 
onator light 0. In addition, such position resolution is 
allowed by the lOnm and better ground state sizes achiev- 
able in magnetic traps on the atom chip 9] . Future work 
will include exact simulations of the near surface poten- 
tial. 



V. SINGLE-ATOM DETECTION 

In the following we will investigate the properties of 
a microdisk, modeled as a high finesse ring resonator as 
described in the previous section, as a single-atom de- 
tector for quantum information processing on an atom 
chip. The scheme we discuss here is based on homodyne 
detection of the phase change of the light at the output 
of the linear waveguide in the forward direction. There 
are several advantages of this scheme over a correspond- 
ing absorption detection, (i) It allows one to drive the 
atom far off resonance, in which case the precise tuning 
of the disk resonator with respect to the atomic transi- 
tion frequency is of minor importance. Cavity tuning will 
be discussed in more detail in Sec. IVII (ii) If the addi- 
tional loss mechanisms discussed in Sec. IIII CI are small 
compared to the disk-waveguide coupling strength, all of 
the pump light will leave the system through the for- 
ward waveguide output. Therefore for most parameter 
regimes a strong signal can be expected, which allows 
the use of standard photodetectors rather than sophisti- 
cated single-photon counters, (iii) The strong output sig- 
nal also provides stability of the detection scheme against 
weak background scattering processes. 

Our detection scheme works as follows. The output 
field of the straight waveguide is mixed with a strong 
local oscillator field at a 50-50 beamsplitter and the light 
intensities in the two beamsplitter outputs are measured 



and integrated over the observation time r to give the 
total number of detected photons Ni and A^2- The signal 
we are interested in is given by the difference jA'^i — A^2|- 
The phase of the local oscillator is adjusted such that this 
difference is zero when no atom is interacting with the 
disk field. The presence of an atom is then inferred from 
a change in this intensity difference. Assuming that the 
detection is shot-noise limited, the signal-to-noise ratio 
S of the atom detection is given by 

S = « 2V7| A„*^o|| sin(</. - Ml (28) 

where (p {(j)o) is the phase of the waveguide output field 
Aout.+ with (without) an atom, Aout,o is the field ampli- 
tude without an atom, and we have assumed that only 
the phase and not the amplitude of the output is changed 
due to the atom. 

The second quantity of interest is the number of pho- 
tons M spontaneously scattered by the atom during the 
interaction time, given by 

M = 2rr/9ii. (29) 

This should be as small as possible in order to mini- 
mize the backaction of the detection process onto the 
atom. Ideally, M ^ 1 corresponds to the limit of non- 
destructive measurements. Finally, we note that S and 
M scale differently with interaction time r. For any given 
parameters, we can thus choose r in such a way to yield 
a certain, fixed signal-to-noise ratio. As an example, we 
will in the following consider the number Mio of spon- 
taneously scattered photons when r is rescaled to give 
S = 10, 

Mio = lOOM/S"^. (30) 

For simplicity we will assume that the cavity is driven 
on resonance with the disk modes, Ac = 0, thereby mini- 
mizing the effect of other, off-resonant modes. The atom 
is assumed to be far off resonance with respect to the 
cavity mode, ^ F, such that the effect of the atom 
is mainly to provide a phase shift of the cavity mode. 
Under these conditions and in the limit of small atomic 
saturation we can derive analytical approximations for 
5, M and Mio El 

S = 4V7|^™|^, (31) 

M = 4r|A,„p^, (32) 

Mio = 25^. (33) 

Note that Mio is independent of the pump power and of 
the atomic detuning. 

In Fig. [SI we show S, M and Mio as a function of the 
input power to the linear waveguide. For a weak pump 
the signal-to-noise ratio increases with power since more 
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FIG. 5: (a) Signal-to-noise ratio versus pump intensity for 
disk diameter 30/im (solid line) and 15/xm (dashed), (b) Cor- 
responding photon scattering M (top curves) and Mio (bot- 
tom). Gap size is 0.3/im, waveguide width is 0.6/im, a = 2nm, 
Lc = lOnm, Aa = lOOF, r = lO/^s, and the atom is assumed 
to sit 50nm away from the disk surface. 
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gap size (|lm) 

FIG. 6: (a) Signal-to-noise ratio S and (b) scattered photons 
Afio versus gap size for weak pumping (|Ai,ip — 10* pho- 
tons/s). Solid curve: D — 3Q/im, dashed: D — 15/im for 
a = 2nm, Lc — lOnm. Dotted curve: D — 15/im, a = Inm, 
Lc = 5nm. Waveguide width is 0.6^m, distance atom-disk is 
50nm, interaction time is 10/is. 



photons are coupled into the cavity and interact with the 
atom. However, because of saturation the atom can only 
interact with a maximum number of photons in a given 
interaction time. Hence, S reaches a maximum value, 
and for stronger pump powers S is decreasing again. The 
number of spontaneously scattered photons M increases 
with pump power and finally saturates at the value Ft. 
Mio shows an approximately linear increase with pump 
power, indicating that the least perturbing atom detec- 
tion for a given value of S is achieved for low atomic 
saturation. There is, however, a trade-off as the required 
interaction time increases in this limit. For low satura- 
tion, pii <ti 1, the numerical results are accurately de- 
scribed by the approximations H31|) - (|33|) . 

Figure shows S and Mio as a function of the gap 
between the microdisk and the waveguide in the limit of 
weak atomic saturation (pn < 0.03 for the shown param- 
eter range). For small enough gap sizes, increasing the 
gap reduces the coupling between disk and waveguide 
modes and therefore the cavity finesse increases. This 
leads to improved signal-to-noise ratios and less sponta- 
neous photon scattering by the atom. For very large gaps, 
on the other hand, the cavity finesse is limited by the ad- 



ditional losses, see Sec. lIII CI In this case, increasing the 
gap even further reduces the number of photons coupled 
back from the cavity into the waveguide and thus reduces 
the detected signal. If the cavity-waveguide coupling ex- 
actly equals the additional losses, no light is transmitted 
through the waveguide at all, which leads to the points 
of 5 = and the corresponding divergence of Mio ob- 
served in the figure. We find numerically that S is maxi- 
mum and Mio is minimum if losses from the disk into the 
waveguide are about 4-5 times higher than the additional 
losses. 

For the parameters of Fig. |^ we find that the minimum 
value of Mio is 0.85 for a disk diameter of D = 30/im 
(solid line) and 0.49 for D — 15/zm (dashed line). The 
reason for this difference is mainly that for the smaller 
disk more energy of the resonant mode is in the evanes- 
cent field. This leads to improved coupling of the atom 
to the mode, that is, to a larger Rabi frequency g as al- 
ready seen in Fig. ^ For both disk diameters, however, 
the minimum value of Mio is below unity, which indicates 
that single atoms can be detected while on average scat- 
tering less than one photon spontaneously. Moreover, we 
observe that this minimum value of Afio is limited by the 
additional losses due to surface roughness. If the surface 
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FIG. 7: Mio versus mode coupling parameter e/nioss for D — 
15/im, gap size 0.5/im, a = 2nm, Lc = lOnm, Aa ~ lOOF, 
l^inp ~ 10* photons/s. Solid curve: Ac — 0, dashed: Ac — e, 
dotted: Ac = — e. 

parameters are decreased by a factor of two to ct = Inm 
and Lc = 5nm, Miq can be as small as 0.13, shown by 
the dotted curve in Fig.EJb). In this case, our atom de- 
tection scheme approaches the limit of a non-destructive 
measurement. 

Note also that the parameters where optimum single- 
atom detection is observed correspond to the strong cou- 
pling regime of cavity QED ^, defined by 5V(kF) > 1. 
In particular, at the minimum points of Miq in Fig. |H| 
we find g^/{Kr) = 48 (L» = iOfim), 73 {D = Ibfim, 
a ~ 2nm, Lc = lOnm), and 290 {D = 15/im, a = Inm, 
Lc — 5nm). 

Finally, in Fig. |7| we show the effect of mode cou- 
pling between the two counter-propagating WGM on 
the single-atom detection. For e smaller than the to- 
tal cavity linewidth k, the main effect of mode coupling 
is to increase the loss rate of the forward propagating 
mode. This decreases the mode quality factor and there- 
fore increases Miq. For e > k, the strong coupling lifts 
the degeneracy between the counter-propagating modes. 
The resulting new frequency eigenmodes of the disk are 
shifted in frequency by ±e with respect to the uncou- 
pled modes. Tuning one of these coupled modes into 
resonance with the pump may thus seem advantageous. 
However, as can be seen from the corresponding curves 
in Fig.[7| this is not the case for some parameter regions. 
Several effects play a role here: (i) our detection scheme 
only measures the forward propagating output of the lin- 
ear waveguide and thus is sensitive to the phase difference 
between the two coupled modes, which in turn changes 
with the frequency shift; (ii) coupling to the backward 
propagating mode can be regarded as an additional loss, 
which shifts the position of the divergences of Mio al- 
ready observed in Fig. |H| (iii) strong coupling between 
the counter-propagating modes fixes their relative phase, 
which creates a standing-wave pattern at the disk surface. 
Thus, the atom detection becomes dependent on the po- 
sition of the atom. In particular, if the atom is trapped 
at a node of the standing wave, no interaction with the 



light occurs and atom detection is impossible. The in- 
terplay between these effects explains the features of Fig. 
[3 Note finally that for high quality disk resonators only 
weak mode coupling is expected, for example, in Ref. [Sj 
a mode coupling parameter e/nioss of 1.5 has been re- 
ported. 



VI. EXPERIMENTAL CONSIDERATIONS AND 
TUNABILITY 

In this section we consider the feasibility of a microdisk 
resonator in a realistic experiment on an atom chip. In 
particular we discuss the effects of finite fabrication tol- 
erances, of temperature fluctuations, and various possi- 
bilities to tune the resonance frequencies of the disk. 

Let us first consider fabrication imperfections due to 
edge roughness, which changes the disk diameter, or due 
to impurities, which change the index of refraction. To 
first order, these uncertainties affect the mode resonance 
through the simple relation fT3| 

Av/v = -An/n- AD/D, (34) 

where Ai', An, and AD are the changes of the resonance 
frequency, the refractive index, and the disk diameter, 
respectively. Assuming typical fabrication tolerances for 
the disk diameter of a few nm and for the refractive index 
of the order of 10~^ therefore leads to frequency shifts of 
the order of tens of GHz. 

Similarly, temperature fluctuations typically give rise 
to changes of the refractive index of the order of 10~^ 
per 1°C. This again will shift resonance frequencies of 
the microdisk by amounts in the GHz regime. Therefore, 
methods to dynamically tune the disk resonator have to 
be implemented on the chip in order to provide frequency 
stability. 

Tunable photonics is also needed if we want to trap, 
manipulate, and measure different kinds of atoms or 
atoms in various internal states. Accurate control of the 
light-matter interaction is especially critical if one is to 
enable the exploitation of quantum technology, for exam- 
ple, in the case of quantum communication in order to 
convert a flying quantum bit (qubit) in the form of a pho- 
ton into a storage qubit (atom), in the case of quantum 
computing where single qubit rotations, e.g., between two 
hyperflne states, are needed, or in the case of sensors to 
measure interference patterns. 

As already noted, high Q devices such as microspheres, 
Fabry-Perot cavities, or microdisks are an extremely ef- 
fective tool for the delicate manipulation and measure- 
ment of subtle quantum states, where a single photon can 
interact many times with the same atom so that signif- 
icant interaction can be achieved. However, such strong 
coupling requires that the device is kept on resonance 
with the exact frequency close to the chosen atom tran- 
sition mill. 

In general, the FSR of a microdisk is orders of magni- 
tude larger than the linewidth, and therefore coincidences 
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FIG. 8: Required change of disk diameter and refractive in- 
dex for a full FSR scan versus disk diameter at A = 780nm. 
AD/D and Anjn are presented by dotted and dash-dotted 
curves for n = 1.454, and by solid and dashed curves for 
n = 2.17, respectively. Typical values for actual materials are 
also presented by the long base arrows. 



between the fundamental transverse WGM and specific 
atomic transitions are extremely unlikely. To keep a 
WGM resonance near the wavelength of interest we thus 
need to change the diameter or the refractive index of the 
disk. In order to achieve any given frequency, the tuning 
range has to be of the order of the FSR. Furthermore, the 
procedure should be stable, reversible, and fast enough 
to compensate for temporal instabilities such as those 
arising from temperature fluctuations of the chip. Under 
these conditions a resonant mode close to the required 
atomic frequency can always be found, even if the disk 
was initially fabricated with some mismatch in diameter 
or in refractive index. 

In order to scan a full FSR we require that /S.v/v « 
where I is the longitudinal mode index. Together with 
Eq. (|34|l this yields the required relative change of disk 
diameter and refractive index. Results for this relative 
change are shown in Fig. |S1 versus the disk diameter for 
fused silica [n — 1.454 at A = 780nm) and for n — 2.17 at 
A — 780nm. Note however that since a few higher order 
radial modes can also be used for some applications, there 
may in fact exist several usable resonances within each 
FSR. 

At least three possibilities may be considered to realize 
a tunable microdisk resonator: UV, electro-optical, and 
piezo tuning. 

UV tuning was reported in [2^. There, it has been 
shown that it is possible to change the refractive index of 
Ge doped silica by up to An = 0.006 using UV radiation. 
The drawback is that this procedure is not reversible. 

In the case of electro-optical tuning, the effect is 
achieved by a uniform electric field that changes the re- 
fractive index and thus allows one to tunc the resonance 
wavelength. The idea is to cover both the bottom and the 
top of the disk with a metal layer, and apply a voltage to 
create the necessary electric field. For example, electro- 
optic crystaline materials usually have a relative change 
of refractive index of 0.01%- 1% for a field of lO^V/m (5V 



for our 5/xm disk thickness). Though demanding a more 
complex fabrication process, such crystals may indeed 
be used. Figure |S1 shows that disk diameters as small as 
15/im will enable a full FSR scan assuming An/n = 1%. 
Other methods for such changes of the refractive index 
also exist. For example, it was reported in 27] that tun- 
ing an InP microdisk can be achieved utilizing free carrier 
injection, which results in An ~ 0.002 up to An ~ 0.01. 

By using a piezo effect one can change the diameter of 
the disk. In this case, one has to fabricate the disk from 
a transparent piezo material. The voltage necessary for 
tuning will be obtained as above by two electrodes evap- 
orated below and above the device. Transparent piezo 
materials (such as BaTiOa) have a piezo electric coeffi- 
cient of order 10~^°m/V which leads to AZJ/D = 0.003 
for electric fields of 30 x lO^V/m (150V for our 5/xm disk 
thickness). Figure |S1 shows that in this case disk diame- 
ters as small as 30/i7n will enable a full FSR scan. For 
both disk size tuning and refractive index tuning, higher 
voltages and more exotic materials should allow for even 
smaller disk sizes. 



VII. SUMMARY AND CONCLUSIONS 

We have investigated the feasibility of initiating strong 
interaction between photons and single atoms on an atom 
chip by utilizing microdisks as high Q resonators. 

Our calculations show that whispering gallery modes 
with high quality factors of up to 10* can be supported by 
disks with diameters of several tens of microns. Single- 
atom detection with high signal-to-noise ratios and prac- 
tically no spontaneous photon scattering can be achieved 
in such a system, and hence non-destructive measure- 
ments may be possible. 

We have also discussed different methods which would 
allow one to tune the resonance frequencies of a microdisk 
over a full free spectral range. 

These results suggest that the regime of strong cou- 
pling between atoms and photons can be achieved in op- 
tical resonators as small as five microns in diameter. This 
could open the way for compact arrays for quantum in- 
formation processing. 
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